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UNSTIFFENED CIRCULAR CYLINDERS SUBJECTED TO 
BENDING AT ELEVATED TEMPERATURES 
By Eldon E* Mathauser and Avraham Berkovits 

SUMMARY 


A method based on a semienipirlcal procedure Is presented for pre- 
dieting static strength and creep buckling of unstiffened circular 
cylinders subjected to pure bending at elevated temperatures. The 
method is applicable to cylinders that are loaded into the inelastic 
stress range prior to buckling and fail in a local mode. The predicted 
bending moments associated with static strength and creep buckling are 
compared with experimental data (tak.en from NACA RM 5TE1T) obtained 
from tests at 500^ F on 5052-0 aluminum-alloy cylinders with radius - 
thickness ratios ranging from 125 to 250 . 


INTRODUCTION 


The development of semiempiricad. methods of analysis applicable to 
the determination of static strength and creep buckling of structural 
con^jonents at elevated temperatures has been the object of several 
studies made by the NASA (formerly NACA). A summary of some of these 
studies is given in reference 1. These studies were concerned with the 
elevated-temperature behavior of various structural components but did 
not include unstiffened circular cylinders. In view of the importance 
of cylinders in structural applications at elevated temperatures^ some 
investigations of cylinder behavior at elevated temperature have been 
made recently by several investigators. Studies of elevated-tenperature 
strength and creep buckling of axially compressed cylinders are reported 
in references 2 and 5- Experimental strength and creep buckling results 
from unstiffened circular cylinders subjected to pure bending at ele- 
vated temperatures are given in reference 

In the present paper, the semi empirical methods of analysis pre- 
sented in reference 1 are applied in the determination of elevated- 
temperat-ure static strength and creep buckling of unstiffened circular 
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cylinders ■ subjected to pure bending and loaied into the inelastic stress 
range prior to buckling. The bending strength of the cylinders is deter- 
mined on the basis of an existing local bucking stress relationship and 
the compressive stress- strain curve for the material. Creep buckling is 
determined in a similar manner with the use of isochronous stress-strain 
curves obtained from compressive creep data. The predicted results are 
compared with the experimental data given in reference 4. The prediction 
of creep buckling stresses from stresses associated with either experi- 
mental or calculated static strength of the cylinders is discussed. 


SYMBOLS 


A area, sq in. 

E Young * s modulus , ksi 

Eg secant modulus, ksi 

E-t, tangent modulus, ksi 

k buckling coefficient 

M bending moment, in -kips 

n total number of elements 

r radius, in. 

t thickness, in. 

y perpendicular distance from axis :>f bending to element of 

area in cylinder, in. (see fig. 3) 

e strain 

a stress, ksi 

local buckling stress, ksi 

T time, hr 

T] plasticity reduction factor 

Subscript : 


i 


pertains to element 
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METHOD OF ANALYSIS 


Static Strength 

In the present study, the determination of the static bending 
strength of cylinders is made on the basis of the following equation 
which defines the elastic local buckling stress for axially compressed 
cylinders : 


J (1) 

where k is the buckling coefficient and has a theoretical value of 
approximately 0-6, E is the elastic or Young’s modulus, and t/r is 
the thickness -radius ratio of the cylinder. Equation (l) was derived 
for axially compressed cylinders (ref. 5) and has often been used to 
determiine local buckling stresses for cylinders subjected to bending. 
Experimental results indicate that the buckling coefficient k is not 
a constant but, rather, is dependent on the radius-thickness ratio of 
the cylinder. (For example, see ref. 6.) Furthermore, for cylinders 
stressed into the inelastic range of the material prior to buckling, 
the elastic modulus E in equation (l) is replaced by an effective 
modulus y]E as follows : 


acr = kT]E J (2) 

Different material parameters have been assumed for the plasticity reduc- 
tion factor T|; for example, E^/E, in reference 2 and Eg/E, in refer- 
ence 5. In the present study, tj is assumed to be given by the param- 
eter ^E^Eg/s, which is a simplified approximation for the plasticity 

reduction factor given in reference 7- Equation (2) may thus be written 
as 


^cr “ k^EtEg i 

where and Eg are evaluated from the material stress -strain curve 

at the stress a 

cr ■ 

Equation (2) is assumed to define the extreme-fiber local buckling 
stress for cylinders subjected to pure bending and is applicable if 
buckling occurs in either the elastic or inelastic stress ranges of the 
material at room or elevated temperatures. The buckling stress is also 
assumed to be equivalent to the maximum extreme-fiber bending stress 





that can be applied on the unstiffened circalar cylinders. This assump- 
tion, in general, is in good agreement with experience. In addition, it 
is assumed that no local deformations occur prior to buckling and that 
the strain distribution is linear over the cross section of the cylinder 
in both the elastic and inelastic stress ranges of the material. 

The solution for the buckling stress of cylinders loaded into the 
inelastic stress range from equation (2) can be accomplished conveniently 
with the aid of a plot such as shown in figure 1- The solid curve deter- 
mined from the material stress -strain curve indicates the variation of a 

with cr/qE, where r\E = yE^Eg. The dashed Line indicates a constant value 

of k The stress defined by the intersection of the solid and dashed 
curves gives the extreme-fiber buckling stress for the cylinder. 

The magnitude of the buckling moment associated with the extreme- 
fiber buckling stress is generally of interest. If the cylinder is 
loaded into the Inelastic stress rajige prior to achieving Ccr^ 
relationship between the applied bending monent and the extreme-fiber 
stress becomes nonlinear as shown, for example, in figure 2. In this 
method of analysis, the magnitude of the monent corresponding to exi 
inelastic extreme-fiber stress is determined by a numerical procedure. 

The numerical procedure discussed in the following paragraph is based 
on the assumption that the tensile and con^ressive stress-strain curves 
are similar. Modifications to this procedure for dissimilar tensile and 
compressive material properties are described subsequently. 

In the numerical analysis, the cylinder is subdivided into kO ele- 
ments of equal area. A quadrant of the cylinder indicating this sub- 
division is shown in figure 5* Forty elements were selected on the basis 
of a preliminary study which indicated that adequate accuracy for the 
purpose of this analysis was obtained with this number of elements. For 
a given extreme-fiber stress, the magnitude of the stresses in all other 
elements is defined from the material stress-strain curve by making use 
of the assumption that a linear strain distribution exists in the cylin- 
der. The magnitude of the bending moment corresponding to the defined 
stress distribution is obtained from the monent equilibrium equation 


cjy dA = M 


( 5 ) 


which is written as follows for use in the lumerical procedure: 


n 

^ cJiYiAi = M 
i=l 


( 4 ) 
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where defines the stress in an element, defines the perpendic- 

ular distance from the axis of bending (see fig. 5) to the element, and 
Aj_ is the area of the element. Equation (4) may then be solved by taking 

C7j_ in the extreme fiber equal to determined from equation (2) . The 

solution for the static bending strength is thus complete inasmuch as fig- 
ure 1 defines the extreme-fiber buckling stress and equation (4) defines 
the buckling moment corresponding to this stress. 

When the tensile and compressive stress-strain curves are not simi- 
lar, an additional equation involving thrust equilibrium is required for 
determination of the relationship between M and a. The equation is 


I 


a 


dA = 0 


(5) 


( 6 ) 


The magnitude of the moment corresponding to a given extreme-fiber com- 
pressive stress is obtained by solving equations (4) and (6) simul- 

taneously by an iterative or successive approximations procedure. This 
method of analysis also defines the neutral or bending axis for the 
cylinder . 


which may be written for the numerical procedure as follows : 

n 

cr^Aj_ = 0 

i=l 


Creep Buckling 

The determination of the bending moments that produce creep buckling 
of circular cylinders is based on the method outlined previously for pre- 
dicting static strength. Isochronous stress-strain ciarves derived from 
material creep data are substituted for the material stress-strain curves 
to predict the extreme-fiber creep buckling stresses from equation (2) . 
The bending moments corresponding to the extreme-fiber creep buckling 
stresses are determined from equation (4) for similar tensile and com- 
pressive creep properties and from equations (4) and (6) for different 
tensile and compressive creep properties. 
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EXPERIMENTAL RESULTS 


Cylinder Data 

The experimental data given in reference k on the elevated-temperature 
static strength and creep buckling of circular cylinders were used as a 
basis for comparison with the results predicted by the methods outlined in 
the present study. The test specimens reported in reference 4 were 
unstiffened circular cylinders 48 inches long and l6 inches in diameter. 
These specimens were fabricated from 5052-0 aluminum-alloy sheets with 
nominal skin thicknesses of 0.064^ O.O5I, 0.040, and O.O52 inch. These 
skin thicknesses correspond to nominal radius-thickness ratios of 125, 

157^ 200, and 25O, respectively. The cylinders were subjected to pure 
bending at room temperature and 500^ F to c.etermine static strength and 
were also subjected to constant moments at 500^ F to determine creep life- 
times. The results of both the strength arid creep tests are presented and 
discussed in the section entitled "Comparis.on Between Predicted and Experi- 
mental Results." 


Material Stress-Strnin Data 

No materials data were given in reference 4 for the 505^-0 aluminum- 
alloy sheets. In order that comparison beiween the experimental and cal- 
culated results could be made, tensile and compressive stress-strain and 
creep tests were performed on this material for sheet thicknesses of 
0.064 inch and 0.052 inch. These tests were made in both the with-grain 
and cross-grain directions of the sheets. The tensile and compressive 
stress-strain results are given in figure ^ . These data were obtained 
from the 0.064-inch-thick sheet tested in ihe with-grain direction. The 
specimens were exposed to the test temperature for 30 minutes prior to 
loading. No significant differences were cbtained from tests in the 
cross-grain direction for this sheet thickress or from tests in either 
grain direction of the 0.032-inch-thick sheet. It was thus concluded 
that the stress-strain properties for the 5052-O aluminum-alloy sheet 
were essentially independent of grain direction and were also independent 
of sheet thickness. 


Material Creep Dsta 

Tensile and compressive creep tests were performed on the 
5052-O aluminum-alloy sheet to provide material creep data for the 
determination of cylinder lifetime. The results of these material 
creep tests are shown in figure 5- Th^ curves in figures 5(a) and 5(b) 
indicate tensile creep results obtained frem sheets tested in the with- 
grain direction at 500° F. Results for 0.C64-inch-thick and 0.052-inch- 
thick sheets are indicated. These material specimens were exposed to 
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the 500^ F test temperature for 50 minutes prior to the application of 
the constant load in a tensile creep testing machine. Similar results 
were obtained from the tensile creep tests on the material in the cross- 
grain direction. 

Figure 5(c) indicates the compressive creep test results which were 
obtained from the 0.064-inch-thick sheet tested in the with-grain direc- 
tion. The compressive creep specimens were tested in a V-groove edge- 
support fixture and were also exposed to the 500° F test temperature for 
50 minutes prior to application of the constant load. In view of the 
similar material properties obtained from the stress-strain tests in the 
with-grain and cross-grain directions, no compressive creep tests were 
performed on the 0.064 -inch-thick sheet in the cross -grain direction or 
on the 0.032-inch-thick sheet in either grain direction. 

A cross plot of the compressive creep data given in figure 5(c) at 
constant times gives the isochronous compressive stress-strain curves 
shown as' solid curves in figure 6. The isochronous compressive stress- 
strain curves indicate the strain produced immediately on application of 
a given stress plus the creep strain obtained at that stress for the 
indicated times. The dashed curve (labeled a-e) shown for comparison 
purposes is the compressive stress-strain curve taken from figure 4 for 
the material at 500^ F. 


COMPARISON BETWEEN PREDICTED AND EXPERIMENTAL RESULTS 


Static Strength 

The predicted extreme-fiber buckling stresses appropriate to the 
four radius-thickness ratios of the test cylinders are indicated in fig- 
ure 7 tbe intersections of the solid and dashed lines. The solid 
curve in this figure indicates the variation of a with a/r]E obtained 
from the compressive stress-strain curve (fig. 4) at 500° F. The dashed 

lines indicate constant values of k i which are applicable to the four 

radius-thickness ratios of the test cylinders. The magnitude of the 
buckling coefficient k was taken from figure 8. This variation of k 
with radius-thickness ratio was evaluated from reference 6 and is appro- 
priate for unstiffened cylinders subjected to bending. 

The magnitude of the bending moments that correspond to the extreme- 
fiber buckling stresses indicated in figure 7 calculated from equa- 
tion (4). These moments are indicated by the end points of the curves 
given in figure 9- The variation of M with a is shown for all values 
of extreme-fiber stress up to The relationship between M and a 
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is linear up to an extreme-fiber stress of approximately 7 ksi, which 
corresponds to the proportional limit stress for the material at 500° F. 
(See fig. 4.) Above 7 ksi, the variation of M with 0 becomes non- 
linear; however, the deviation from linearity is small, particularly for 
the cylinders with the high radius -thickness ratios. 

A comparison between the experimental and predicted strengths of the 
cylinders at 500° F is given in figure 10. For radius -thickness ratios 
of 157 and 200, satisfactory agreement between the experimental and cal- 
culated results was generally obtained. Pc_)r radius-thickness ratios of 
125 and 250 , the predicted moments exceed uhe experimental values. No 
attempt was made to compare the experimental and predicted strengths at 
room temperature. The room-temperature stress-strain curve is character- 
ized by a flat top at the initiation of inelastic action in the material. 
For this type of stress-strain curve, the magnitude of the inelastic 
extreme-fiber buckling stress is defined; however, the magnitude of the 
corresponding strain is indeterminate. In view of the indeterminate 
magnitude of the strain in the extreme fiber, the stress distribution 
over the cross section of the cylinder is unknown and the use of equa- 
tion (4) for the determination of the bending moment is precluded. 


Creep Buckling 

The bending moments which produce creep buckling of the cylinders 
for radius-thickness ratios appropriate to the test data of reference 4 
were predicted from the method outlined previously for determination of 
static strength. Isochronous stress-strain curves were substituted for 
the material stress-strain curves in the ce Iculations . The initial cal- 
culation was based on compressive creep daia only amd is described in the 
following paragraph. Calculations based or tensile creep data only as well 
as calculations based on both tensile and compressive creep data were also 
made and are discussed subsequently. 

The extreme-fiber compressive stresses that produce creep buckling 
i-i various times were obtained from figure 11. The solid curves indicate 
the variation of 0 with 0 /t]E obtained from the isochronous stress- 
strain curves shown in figure 6. The dashed lines in figure 11 are con- 
stant values of k i which are applicable to the indicated radius- 

thickness ratios. (The values of the buckling coefficient k were 
taken from fig. 8.) The bending moments corresponding to the extreme- 
fiber creep buckling stresses are indicated by the end points of the 
curves in figure 12. The bending moments were calculated from equa- 
tion (4) and the isochronous stress-strain curves in figure 6. The var- 
iation of M with 0 is shown for all values of extreme-fiber stress 
up to 0cr* Although the relationship between M and 0 is nonlinear. 
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the deviation from linearity is small and compares approximately with 
the deviations obtained from the static strength calculations shown in 
figure 9* The creep buckling moments indicated in figure 12 are com- 
pgired with the experimental results in figure IJ. For lifetimes less 
than 0.2 hour^ the experimental creep buckling moments are generally 
less than the predicted moments. For longer lifetimes^ the experimental 
moments exceed the predicted moments. 

Calculations for creep buckling moments based on tensile creep data 
only were made for purposes of comparison with the previous calculations 
based solely on compressive creep data. The buckling moments predicted 
from the tensile creep data parallel the results obtained from the com- 
pressive creep data shown in figure 15 ^ and^ in general, were approxi- 
mately 20 percent less than the buckling moments obtained from the com- 
pressive creep data alone. These results were anticipated on the basis 
of a comparison between the tensile and compressive creep data plotted in 
the form of a against a/rjE as shown in figure 11. The tensile creep 
data yielded values of predicted buckling stress that were approximately 
20 percent less than those obtained from the compressive creep data. 

Calculations of creep buckling moments made on the basis of both 
tensile and compressive creep data were also made with the use of equa- 
tions (2), ( 4 ), and ( 6 ). These predicted creep buckling moments again 
parallel the results obtained from the compressive creep data alone. The 
predicted buckling moments in this case were approximately 5 percent less 
than the corresponding buckling moments shown in figure IJ. These calcu- 
lations based on both tensile and compressive creep data did not improve 
the correlation between experimental and calculated creep buckling moments 
over that obtained from compressive creep data only. Furthermore, the 
computation for creep buckling was simplified considerably when made on 
the basis of compressive creep data only. Analysis for creep buckling 
based on both tensile and compressive creep data may be necessary only 
when considerably different tensile and compressive creep properties 
exl st . 


CONCLUDING REMARKS 


An analysis based on a semiempirical procedure has been presented 
for predicting the elevated-temperature static buckling strength and the 
creep buckling strength of unstiffened circular cylinders subjected to 
bending. The determination of buckling stresses by this procedure is 
based on the observation that a shell structure of a given geometry 
appears to buckle at a stress corresponding to a constant value of a/r|E 
regardless of material or test tenperature, provided that a suitable 
expression for the plasticity reduction factor r\ is chosen (a = Stress; 

E = Yo\mg*s modulus). Comparison of the predicted results with experi- 
mental results (taken from NACA RM 57E17) obtained on 5052-0 aluminum-alloy 
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cylinders with radius -thickness ratios ranging from 125 to 25 O indicates 
that this procedure is valid for determining the buckling stresses asso- 
ciated with either static strength or creep buckling of cylinders. Iso- 
chronous stress-strain curves are employed in this procedure for deter- 
mination of the creep buckling stresses. 

In view of the satisfactory correlation between the predicted and 
experimental static strength and creep buckling results^ it appears that 
creep buckling stresses may be predictable solely from a knowledge of 
the static strength of the cylinders at the same or lower test tempera- 
tures and the appropriate material stress-strain and material creep data. 
Additional tests of cylinders fabricated from other materials and tested 
at several temperatures will help to determine the validity of this 
assumption . 

Langley Research Center^ 

National Aeronautics and Space Administration, 

Langley Field, Va., March 24, 1959- 
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Figuz*0 1 .— Gr 3 .phica.l 501111,1011 lor local buckling slrcss. 



Figure 2 .- Example of variation of appliec. moment with extreme-fiber 
stress for cylinder loaded into inelastic stress range prior to 
buckling. 
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ore 4.- Tensile and compressive stress-strain curves for 5052-0 
aluminum-alloy sheet. Thickness* 0.064 inch; nominal strain rate 
0.002 per minute; exposed at 500^ F for 50 minutes prior to load! 
in 500^ F tests; tests performed in vith-grain direction of sheet 
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(a) Tensile creep curves; sheet thickness, 0-064 inch. 



(b) Tensile creep cuxves; sheet thickness, 0.0J2 inch. 



(c) Compressive creep curves; sheet thickness, 0.064 inch. 

Figure 5.- Material creep curves for 505^-0 aluminum-alloy sheet at 
500 ^ F. Exposed at 500^ F for 50 minrtes prior to loading; tests 
performed in with-grain direction of sheet. 



Figure 6 .- Isochronous compressive stress-strain curves for 
5052-0 aluminum-alloy sheet at 500 ^ F. 



Figure 7 *- Graphical solution for extreme-fiber buckling stresses for 
5052-O aluminum -alloy cylinders at 500 ^ F. 



Figure 8*- Variation of local buckling coefficient with radius- 
thickness ratio for unstiffened circi:lar cylinders subjected 
to bending. (Evaluated from ref. 6.) 
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Figure Variation of applied moment witi extreme-fiber compressive 
stress for 5052-0 aluminum-alloy cylinders at 500^ F. 




Figiire 10.- Comparison between calculated and experimental buckling 
moments for 5052-0 aluminum-alloy cylinders at 500° F. 



Figure 11.- Graphical solution for extreme-fiber corr 5 )ressive stresses 
that produce creep buckling of 5052-0 aluminum-alloy cylinders at 

T? 


Figure 12.- Variation of applied moment with extreme-fiber compressive 
stress for 5052 -O aluminum-alloy cylinders at 500° F. 



Figure ij.- Comparison between predicted and experimental creep buckling 
moments for 5 O 52 -O aluminum-alloy cylinders at 500° F. 
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